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We study the Laplacian operator of an uncorrelated random network and, as an application, 
consider hopping processes (difltusion, random walks, signal propagation, etc.) on networks. We 
' develop a strict approach to these problems. We derive an exact closed set of integral equations, 

' which provide the averages of the Laplacian operator's resolvent. This enables us to describe the 

, propagation of a signal and random walks on the network. We show that the determining parameter 

OA ' in this problem is the minimum degree qm of vertices in the network and that the high-degree part 

f— I ' of the degree distribution is not that essential. The position of the lower edge of the Laplacian 

, spectrum Ac appears to be the same as in the regular Bethe lattice with the coordination number 

Qm- Namely, Ac > if gm > 2, and Ac = if gm < 2. In both these cases the density of eigenvalues 
■ p (A) ^ as A Ac + 0, but the limiting behaviors near Ac are very different. In terms of a distance 

from a starting vertex, the hopping propagator is a steady moving Gaussian, broadening with time. 
This picture qualitatively coincides with that for a regular Bethe lattice. Our analytical results 
f~i ' include the spectral density p(A) near Ac and the long-time asymptotics of the autocorrelator and 

O ' the propagator. 
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INTRODUCTION 



I , The Laplacian spectra of random networks determine a wide circle of processes taking place on these networks, 
' see, e.g., [l|, 0, S S IE S S H, and references therein. Random walks, signal propagation, synchronization, and 
^ [ many others are among these processes. This is why the problem of Laplacian spectra of random networks (especially, 
O , of its low-eigenvalue part which determines the long-time behavior of relevant processes) is considered as one of 
central problems of graph theory and the science of complex networks. In this paper we essentially resolve this 
problem applying the strict statistical mechanics approach to uncorrelated random networks with arbitrary degree 
distributions. These random graphs constitute a basic class of complex networks. 
^ [ One should note that leading contributions to the spectra and the asymptotics of the random walk autocorrelator 
were found by Bray and Rodgers in 1988 in the particular case of the Erdos-Renyi graphs 3]. It is important that 
these classical graphs necessarily have dead ends and vertices with two connections. We will show that the absence 
I I of these vertices in a network qualitatively changes the spectra and the random walk asymptotics. Random walks 
• ' on hierarchically organized, deterministic, scale- free graphs were studied by Noh and Rieger in Ref. p^ . Due to a 
, very specific organization of these graphs, their results are not applicable to equilibrium networks. This is also the 

' case in respect of the recent numerical work of Kujawski, Tadic and Rodgers who found the autocorrelator of 

, a random walk on a growing scale-free network by performing extensive numerical simulations. Their network was 
I!" ■ strongly correlated in contrast to the configuration model of a random graph, which we use in this work. 
. ^ ] For the sake of clearness, let us remind basic notions and terms for random networks. For more detail see [ll|, [l^, 
■ [H, [13, [iM EB, [i3|- ^ graph is completely defined by its x adjacency matrix A, whose elements Aij are the 
^ I numbers of edges between i and j. The vertex degree of vertex i is the number of edges, attached to this vertex: 
- - ' qi = X^jLi -^ij — -^ji- random networks, Qi is a random variable with a degree distribution H (q) = {S {q — qi)). 

In traditional mathematical models, Ii{q) is a rapidly decaying function with a well-defined scale. For example, in 
the Erdos-Renyi model 18], which is a standard one, Il{q) is a Poisson distribution decaying as [q/eqY , i.e., faster 
than any exponent. In contrast to these models, in most of real-world networks degree distributions are heavy tailed. 
After the work [T9| . they are usually approximated by a power-law ^ q^^ in the range of sufficiently high degrees. 
Note that the validity of this fitting is limited because real-world networks are small (even the WWW has only about 
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TABLE I: Asymptotics of the Laplacian spectral density p(A), autocorrelator Pq (t) and propagator Pi (t) for the random 
uncorrelated networks where 11 (q-m) is essentially distinct from and 1. Here P;*^'' = Piit ^ oo) are stationary values of the 
correlator given by Eq. (I40p for I = and Eq. (|47p otherwise. P — n {qm — 1)^''* In {qm — 1). The values of the parameters in 
the pre-exponential factors are ^ = 9/10 and 4/3, rj = —7/30 and 1/18, and — 13/30 and 5/18 for qm ~ 1 and 2, respectively. 
V and D are determined by the full form of the degree distribution 11 (q). 





Minimum vertex degree qm > 2 


Minimum vertex degree Qm ~ 1 or 2 


Spectral edge Ac 


qm - 'i'^/qm — 1 





Spectral density 


«^P[2^ dexp(^)],Eq. 113 


p'l^'^^S (A) + const A-« exp (^-a/V\j , Eqs. IMJ, IZll) 


Autocorrelator 


exp [-Act - ff-tj In^t] , Eq. (|58|l 


p^,"''' + const t" exp [-3 (f )'''^ t'''''] , Eqs. ^ 


Propagator at Z ~ t 


^exp[-iiz#],Eq. m 


^ exp [-i^] , Eqs. §1, m 


Propagator at t^oo 


fiQ (-Ac) exp [-Act - P^t/ In^t] 


pI"""^ + ct"'' exp [-3 (t - Z/i;)'/='] , Eqs. (IH2J 



10"'^'^ vertices), and so high degrees are not observable. It is commonly beUeved that the "scale- free networks" are 
greatly distinguished from the others in every aspect. This widespread behef actually implies a division of all networks 
into two classes: "scale-free networks" and all others. In contrast to these beliefs, we here show that scale-free (or, 
more generally, heavy tailed) architectures of networks are not essential for a lower edge of the Laplacian spectra and 
the long-time behavior of random walks characteristics. The resulting dependences are determined by the minimum 
degree of vertices in a network. Heavy tails determine some coefficients and amplitudes but not a type of these 
singularities. 

In this paper we study properties of the Laplacian operator 

Lij — QiSij — Aij (1) 

on an uncorrelated random network near the lower edge of its spectrum, and, respectively, the hopping motion of 
some carrier ("signal") from one vertex to another at large times. This operator corresponds to the process described 
the following dynamic equations for the probability pij {t) that at time t a particle is at vertex i if at time it was at 
vertex j : 

N 

Pij (*) ^ X! ^^kPkj{t) - QiPijit), PtjiO) = dij. (2) 

This is a random walk where the rate of hopping along any edge is set to one. Other versions of the Laplace operator 
and corresponding processes, which are also widely discussed in literature, are listed in Appendix [K\ 

We use the configuration model of an uncorrelated network [20l [2l| , which is a maximally random network with 
a given degree distribution. It is convenient that (i) this model is statistically homogeneous, (ii) all its vertices are 
statistically independent, and (iii) it has a locally tree-like structure. We consider only infinite networks, that is, 
first we tend the total number of vertices N to infinity {the thermodynamic limit) and only afterwards study network 
characteristics. If, say, we study a random walk, then a particle should be still much closer to an initial vertex than 
the diameter of the network ~ IniV. In other words, we consider the process at so short times that the number of 
vertices, where the walking particle may be found, is negligible compared with the network's size N. We will see that 
this imposes strong limitations to the applicability of our results due to the "small world" feature of the networks 
under consideration. 

We will show that for the Laplacian spectra and for random walks, the crucial property of the random uncorrelated 
network is the minimum degree qm of its vertices. We suppose that the value of the degree distribution at qm essentially 
differs from and 1. We also assume that Qm > 0, because the contribution of isolated vertices is trivial. Our results 
are summarized in Table Hand in Fig. [T] Note an unusual singularity of the spectral density in the case > 2. 

As is natural, the calculation of the spectrum is reduced to the study of the trace of the Laplace operator's resolvent. 
To describe the propagation of the signal in the network, one must know the non-diagonal elements of the resolvent. 
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FIG. 1: Laplacian spectral density p (A) for networks with different minimum vertex degree qm'- (a) g™, > 2, (b) qm = 2, and 
q-m = ^ (the vertical line at A = represents a (5-function peak). 



Here we calculate the asymptotics of their average values. It allows us to obtain the time and distance dependences 
of the signal's propagator pij (t) — pi (t) when the distance between initial i and final j vertices, / is much smaller 
than the diameter of the network, I ^ InN. 

Why is the minimum vertex degree so important in these problems? Note that in respect of random walks and 
Laplacian operator related problems, infinite uncorrelated networks are equivalent to infinite Bethe lattices with 
coinciding degree distributions. (Recall that a Bethe lattice is an infinite tree without borders.) Let us compare 
two Bethe lattices — random, with the minimum coordination number qm, and regular, with the coordination number 
equal to Qm- It is clear that the autocorrelator in the random Bethe lattice cannot decay slower than in the regular 
Bethe lattice with this Qm- If Qm > 2, then in this regular Bethe lattice, pa {t) = po (t) ^ t~^/'^ exp (— Aci), where 

Ac = qni - '^y/q-in - 1- (3) 

Ac is also the spectral boundary in the Laplacian eigenvalue density p{X) of this regular Bethe lattice, where p{X) ^ 
VA — Ac, near Ac. Thus, the spectral boundary for an infinite uncorrelated network in principle cannot be lower 
than that for the regular Bethe lattice with the same qm- Moreover, these borders coincide. The reason for this is 
the following feature of the configuration model of an uncorrelated network. Let the number of vertices N in this 
model approach infinity. Then the mean number of given finite regular subgraphs with coordination number qm 
grows proportionally to N. We stress that although this number rapidly decreases with a size of these subgraphs, 
it is proportional to N for any given subgraph size. In the arbitrarily large subgraphs, the lowest eigenvalues are 
arbitrarily close to the spectral boundary of the corresponding regular Bethe lattice. The number of these eigenvalues 
is proportional to the number of these subgraphs and so proportional to N . Now recall that the total number of 
eigenvalues in the spectrum is N . Therefore, indeed, the spectral borders for the configuration model and for the 
regular Bethe lattice with qm coincide. 

The statistics of these regular tree subgraphs determine the singularity of the resulting p{\) at the edge Ac- The 
rapid decrease of of the number of these subgraphs with their size results in specific singularities, with all derivatives 
zero, represented in Table [D 

The random networks with qm — 1,2 markedly differ from those with g,„ > 2. In the configuration model with 
qm — 1, 2, chains and chain-like subgraphs are statistically essential. Let us first discuss the case qm = 2. The Bethe 
lattice with coordination number 2 is a usual infinite chain. It has the spectral boundary Ac = 0. Near this edge, 
/9(A) ^ A~^/^. Thus, the edge of the spectrum of the uncorrelated network with = 2 is zero. We will show that 
the statistics of chain subgraphs in this configuration model differ from those for the case > 2. This results in 
different asymptotics presented in Tableland, schematically, in Fig. [T] 

If = 1, chain- like subgraphs are also present in the configuration model. These are, however, more chains (see 
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FIG. 2: Chain with finite tree-hke branches in a random network with the minimum vertex degree 



= 1. 



Figl2]) with branches attached. Nonetheless, these subgraphs result in the spectrum edge Ac = and in the same 
asymptotics as for = 2. When = 1, numerous finite components components are present in the network. Their 
mean number is proportional to N. Each of connected components gives one zero eigenvalue in the spectrum. This 
leads to a J-function peak at A = in the spectral density. 

The found singularities of p(A), with all derivatives zero, have a direct consequence for observations in finite 
networks. Even in a huge uncorrelated network, the observed minimum eigenvalue A2 will be far from the spectral 
edge Ac predicted for an infinite network. Let us roughly estimate A2(A^) based on the spectral densities from Table [J 
The condition N J^^^^"* dXp{X) ^ 1 leads to the following dependences in the range of large N. If q„i = 1, 2, then 



Thus the approach of A2 (N) to Ac is extremely slow. Note that a very slow convergence of A2 was recently observed in 
the numerical work of Kim and Motter, Ref . [22l | , in which A2 and Ac were compared for networks up to 4 000 vertices. 

In Sec. HI] we strictly formulate the problem. In Sec. IIIII we derive a basic set of integral equations. Solving these 
equations enables us to obtain the Laplacian spectrum p (A) for uncorrelated random networks and to describe the 
random walk on the networks in the thermodynamic limit. In Sec. IIVI we study the final value of the propagator 
^(cq) — _> 00), which is the equilibrium probability to find a signal at distance I from a starting vertex. We 

describe p^^'^^ in terms of I and of the degree distribution 11 (q). Furthermore, we find the coefficient of the S (A) term. 
In Sec. |V] we present general solutions of the integral equations of Sec. IIIII and analyse them in three distinct cases: 
qm > 2, q„i = 2, and qm — 1- In Sec. IVII we summarize our results and methods and discuss conditions for their 
applicability. Technical details are given in Appendices. 



The problem of the Laplacian spectrum of a random network is completely equivalent to that of the time dependence 
of the averaged autocorrelator Pq (t) = {Pa (t)) for a random walk. This autocorrelator is the probability that a 
particle returns to the starting vertex after a time t. This quantity is related to the eigenvalue density 



(4) 



and if > 2, then 



A2(A^) - Ac - (lnlniV)-2. 



(5) 



II. FORMULATION OF THE PROBLEM 




(6) 



in the following way: 




(7) 



where Xk are (nonnegative) eigenvalues of the Laplace operator on the network: 




(8) 




(9) 
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We assume, that a particle moves from vertex to vertex by hopping along edges. To every edge we ascribe a hopping 
rate Wij, which is the probability to move from vertex j to vertex i per unit time. Hopping rates are assumed to be 
symmetric and equal 1 for every edge, Wij — Wji — Aij = or 1. In this paper we fix Wij but not the escape rate of 
a particle from a vertex, see Appendix |^ where other forms of a Laplace operator are listed. It turns out that our 
main conclusions are also valid if the escape rate from a vertex is fixed. This case will be discussed in detail in our 
next works. Assume that at t = the particle is at vertex j. Its motion is governed by the master equation for the 
propagator, which is the probability pij (t) that at time t the particle is at vertex i, 

N N N 

Pij (*) = X! ["^'fcPfei (*) " WkiPtj {t)] = ^ Aij [pij (t) - Pj, {t)] = ^ AikPk] (t) - QiPij {t) . (10) 
k=l j=l k=l 

This equation is supplied with the initial condition pij (0) = 6ij. What is the value of the probability 



that at time t the particle is at distance d{i,j) = n from a starting vertex? (The distance is the minimum shortest 
path between two vertices.) Here (• • • ) means the average over some statistical ensemble of graphs (over that of the 
configuration model in our case). 
In the Laplace representation. 



P,fc(s)= / dtp,k{t)e~'\ (12) 
Jo 

the propagator is the resolvent of the Laplace operator: 

P{s)=(s + Ly\ (13) 



Consequently, the density of eigenvalues is expressed in terms of the analytic continuations of the averaged values of 
the autocorrelator: 

pW^T^ [Po (-A - lO) - Po (-A + ^0)] . (14) 
2m 



The inverse relation is 



/+ioo+5 , _ /.oo 

— e^*Po(s)- / dXe'^'piX). (15) 



III. MAIN EQUATIONS 

We assume the thermodynamic limit: N — > cx), and the fraction of vertices with a degree q, N (q) /N — > H(g). 
Here H (q) is a given degree distribution with a finite second moment, q^Il (q) < oo. In this limit, almost all finite 
subgraphs are trees, i.e., they have no closed loops within. The network is uncorrelated, i.e., degrees of any pair of 
vertices, connected or not, are independently distributed random variables. These features allowed us to describe the 
statistics of intertvertex distances [2j|. The problem under consideration is actually related to that work. 

The equation for the resolvent of the Laplace operator p3)) is 

sPtk = 5jfc + ^ Aij {Pjk - Pik) = Sik + ^ A^jPJk - qiPik- (16) 
3 j 

Without lack of generality we choose the initial vertex j ~ 0. 

By definition, the n-th connected component of a vertex i is a subgraph, containing all vertices j within the distance 
d(i,j) < n from the vertex i. For any finite n, in an infinite graph almost any n-th connected component of vertex 
is a tree. Actually, we analyse a random Bethe lattice. Degrees of its vertices are independent random variables. 
Its arbitrary chosen central vertex has the vertex distribution function H {q) — {5 {qq — q)). The other vertices have 
degree distributions equal among themselves but different from H (q) . Non-central vertex i of a degree qi has one edge 
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X/ \l_ 



1,1 nj 




FIG. 3: Vertex {n+l,i) at distance n from the starting vertex 0, its "ancestor" {n,j), and its "descendants" {n + 2,k). 
k = 1,2, . . . q„+i,k — 1, <ln+i,k is the degree of the vertex (n + 1, i). 

directed to the central vertex and hi = qi — I > edges directed from it. Here b is the branching number of the vertex. 
Its distribution is given by 

N 



ni (^) = ^ E (^^^^ - ^)) = (& + 1) , (17) 



where q — 2L/N. 



N N 



^E^'.^^E* (18) 



2 ^ ■' 2 

is the total number of edges in the graph. In Eq. ([TT]) we assume that N oo, L oo, but 2L/N q — J2q 1^ 
where q is some finite number. Hi (6) is the probability that a randomly chosen end vertex of a randomly chosen edge 
in the graph has b = q — 1 edges apart from the chosen edge itself. It is convenient to use distributions 11 and Hi in 
Z-representation (see Appendix IB|) . 

Let {n,i) and (n + be two vertices connected by an edge and at the distances n — 1 and n, respectively, from 
the starting vertex 0. We introduce the following random variable (see Fig. [3]): 



Pn,t j s) - Pn+ l,] (s) 



Tn,, is) = -'^ .T"^ ■ (19) 



It is obvious that the statistical properties of this variable are independent of the particular choice of vertex i in 
the n-th shell of the initial vertex 0. The graph ensemble under consideration is completely defined by the degree 
distribution. All graphs with a given degree distribution have the same statistical weights. This, in particular, implies 
the statistical homogeneity of the ensemble. First, we randomly choose vertex 0. Second, we label all other vertices 
by two indices: the first one is the distance from vertex (the shell's number), and the second index labels vertices 
within the shell. Third, we consider P„i, which is the matrix element of the resolvent for the pair — vertex and vertex 
i at distance n from vertex 0. It is a fluctuating random variable but its statistical properties are independent of the 
choice of i, because every averaging includes averaging over all vertices in shell n. The other fluctuating quantity in 
Eq. (jl9p . P„_|_ij , is, of course, correlated with Pn,i. Nonetheless, due to the statistical independence of the vertices, 
this correlation is independent of the particular choice of the connected pair of vertices. Therefore, we can define the 
distribution function of Tn^ij, which is independent of In the Laplace representation this distribution is defined as 

r„ (s, x) = (exp [-XTn^ij (s)]) . (20) 

Now let us recall that in the infinite network all finite connected components are trees. Moreover, in the thermo- 
dynamic limit the statistical properties of all r variables are the same, i.e., they are independent of n too (see more 
detailed discussion in Appendix[C|. It implies the following important consequences, (i) In the thermodynamic limit, 
i.e., for an infinite network, T„ = T is independent of n. (ii) It is possible to obtain the closed equation for T (s,x). 
(iii) The density of eigenvalues p (A), and, consequently, the autocorrelator po (t) can be expressed in terms of T {s, x) 
(see Appendix [Cjl . 

Equation for T {s,x) may be written as 

e-T {s,x)^l + V^ r (2Viy) e-(i+^)>i [T {s, y)] , (21) 
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where Ii is a modified Bessel function, and (pi (z) is tlie degree distribution of non-central vertices branching numbers 
in Z-representation (see Appendix [B]) . As Res > 0, this function has a solution with all properties of the Laplace 
transform of the distribution density of a nonnegative random variable. This statement may be proved by using an 
approach of Ref . [13] . 

The function ifi{z) has the following properties: (i) ipi (0) = Hi (1), which is the concentration of vertices with 
degree 1 ("dead ends"), and (ii) if the degree distribution Hi (q) decays faster than any exponent for q ^ oo, then 
z = 1 is a point of singularity of ipi (z). The function ipi (z) in the complex plane is analytic within the circle \z\ < 1, 
ifi (1) = 1. The parameter ipi (0) is the crucial one in the division of the graph into connected components, see 
Appendix [Bl We show in this Appendix that the autocorrelator Pq (s) in the Laplace representation is given by 

/>oo 

Po{s) = / dx e-'-v[T{s,x)], (22) 
Jo 

where ip{z) is the Z-transformation of the degree distribution Tl{q). The functions ip and ipi are connected as tpi (z) = 
tp' (z) /(f' (1), so that (fi (1) — (fi (1) — 1. The density of eigenvalues p (A) and time-dependent autocorrelator pQ (t) 
can be obtained from Eqs. (fT4|) and (fT5|) . respectively. 

The propagator p„ (t) at n > may be expressed in terms of some functions C/„ (s, x), for which we have a linear 
recursion, relating J7„ to C/„_i (see Appendix |D]) . These functions are introduced in the following way. Let us choose 
two vertices (n, i) and {n + 1, j), connected by an edge (Fig.[3l). We define 31^}^^ (s) as 

Here the summation is over all those vertices at a distance n + 1 from vertex 0, whose shortest path to vertex runs 
along the edge {n,i) ^ {n + l,j). In other words, the sum in Eq. is over all vertices of the ^-th generation of 
the branch beginning from a chosen edge. For example, sl^\ j (s) = X]fc=Y^ Pn+i,k (s), as one can see from Fig. |31 

Due to the statistical homogeneity of the network ensemble the statistical properties of random variables 5'^''!^ ^ are 
independent of the choice of vertices i and j, if they are connected by an edge. For an infinite network, this statistics 
is also independent of n. The recursion relation can be derived for the following averaged quantity which depends 
only on /: 

Ui {s,x)^ (Sl^l^ is) exp i-xr^,,, {s)] ) . (24) 
The recursion relation is derived in Appendix [D1 It is of the foUowing form: 

Ui (s, x) = / dy lo (2^) e-(i+^) V'l [T {s, y)] Ui-i {s, y) , (25) 







where /q is a modified Bessel function of zero order. This recursive relation is supplied with the initial condition: 

Ui{s,x)^(i + ^\t{s,x). (26) 



dx ^ 

Finally, the Laplace-transformed propagator Pi (s) is expressed as 

/•oo 

Pi{s)= / dxe-''-'p>'[T{s,x)]Ui{s,x). (27) 
Jo 

Equation (|25p may be presented in the form: 

Un^MUn-l, (28) 

where M is a hnear integral operator. Let (s) and tpm (^7^) be its eigenvalues and eigenfunctions, respectively: 

poo 

fi^ (s) e-^ / dy lo (2^) e^fi+^^V'i [T {s, y)] (s, y) = («, x) . (29) 







Note that the general theory of integral operators is usually formulated in terms of which are called their 

eigenvalues. Operator M becomes Hermitian after the substitution tpm (x) — e'^^^^ {tp'i [T (s, a;)]}^^^^ x (*; 2;). The 
kernel of the integral operator in Eq. (|^^ is bounded '2^ if 

/•oo 

dx / dyl^ {2^) e-(2+^)(x+i.)^. {s, x)] p>\ [T (s, y)] < 00. (30) 
Jo 
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This condition is always satisfied if s > 0. Therefore, according to theorems about integral equations with a Hermitian 
bounded kernel [1^, all eigenvalues of this operator are real and finitely degenerate. They form a discrete sequence 
bounded from below, without any condensation point except /i = co. Eigenfunctions are orthogonal and normalizable 
with the weight function e~^^Lp\ [T (s,2:)]: 



dx Lp\\T {s,x)\e '"^ljjkis,x)lpm{s,x) = 6kfn- (31) 

Hence the solution of the recursive relation may be presented as a series in the complete orthonormal set {4>m}, 

Ul (S, X)=Y,A^ (S) (S) (S, x) . (32) 



Taking into account the initial condition (|26p and the orthonormality condition pT|) . the coefficients in this series 
may be written as 

A„,(s)=^ dxi:m{x)^[[T{s,x)]e-'^ [l^^T{s,x). (33) 

Substituting Eq. ((32|) into Eq. (|27p . we get for / > the following relation: 

Pi(s)=^A™(s)B™(s)/i,i„7'(s), (34) 

m 

where 

B„,{s)= / dx e-'^^'\T{s,x)]^„,{s,x). (35) 
Jo 

The resulting propagator P; (s) satisfies the condition of the conservation of the number of particles/signals, which 
in the Laplace representation is Yl^o-^i i^) ~ ^1^- Taking into account Eq. (|34p gives the following form of this 
condition: 

P ( X , A«(s)Bm(s) 1 ,„„^ 



IV. CONTRIBUTION OF FINITE CONNECTED COMPONENTS 



When the minimum vertex degree in the uncorrelated network qm > 2, then (in the thermodynamic limit) the 
network consists of one connected component. If, however, q,n = 1, i.e., 11(1) = IIi (1) = ip' (0) = qipi (0) 7^ 0, then 
connected components exist even in the thermodynamic limit. Their contribution to the propagator at < — > 00 is 
obvious, and can be calculated in a straightforward way. We, however, find this contribution by using the technique 
described in Sec. Illll for the sake of illustration. 

We set in T (s, x) the limit s and x 00, with sx fixed, assuming that there exists a limiting function 

e(z) = lim T{s,z/s). (37) 

s — >-oo 

In Appendix IE 31 we derive the fohowing equation for B: 

e (z) = e- Vi [e (z)] . (38) 

Comparing this equation with Eq. (|B6p from Appendix [Bl one can conclude that Q (z) = H (e^^), that is ip [Q (z)] — 
ip [H (e~^)] = (exp {—zMi)). Here Mi is the size of the connected component with a randomly chosen vertex i. It is 
obvious that the giant connected component, whose size is ~ N , does not contribute to O (2) at any z > as iV ^ 00. 
From Eqs. and ([57)) . we obtain a clear result for the limiting value of the autocorrelator, 

p^"^'^ ^ po (t - 00) = Ihn sPo {s) = dz^[Q{z)\^l^-^y (39) 
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It means that the equihbrium distribution of the signal is homogeneous within its connected component. Passing 
from the variable z to x = (z) and using Eq. (j38p . we calculate this integral: 



j^^'^^ = ^(i,)-iW(tc). (40) 

This result also has a different meaning, 



\i=l / \ Clusters / 

where Nc is the total number of finite connected components. 

A nonzero equilibrium value of the autocorrelator indicates that the degeneracy of the Laplacian eigenvalue A = 
is ~ N. The eigenvectors of this eigenvalue may be chosen in the following way. Each such eigenvector has unit 
vector components in one connected component and zeros in all others. The degeneracy is equal to the total number 
of connected component in the network. 

The t — > oo contribution of finite connected components to pi (t) at Z > may be extracted from the functions: 

ui{z) = lim Uiis,z/s). (42) 

s — *-oo 

In this limit the recurrent relation (j25p turns into 

ui (z) = e-^i [0 (2)] ui-i (z) (43) 
(see derivation in Appendix IE 3|) . From Eq. ([26|) it also follows that ui (z) = Q (z), so 

u,iz)^{e-^^[[e{z)]y-'eiz). (44) 

Let us calculate — Pi{t — oo). The stationary value of (t) at t — > cx) is equal to l/Mj, where Mj is 

connected component with an initial vertex j. Consequently, 

where Q^^^ is the number of vertices at distance / from vertex j. Using Eqs. (f27|) . (|42j) and (|44|). we get 

poo 

pl^q) ^ dx e-^e {x) (fi' [e (x)] {e-^^'i [e (x)]}'"' . (46) 
Jo 

At large I, the region z <C 1, where 8 (z) is close to t^, gives the main contribution to the integral in Eq. (|46p . As a 
result, at large I, we have 



Here we used that Q (0) = tc = (pi {tc) and —8' (0) = tc/ [1 — f'l (tc)], which follows from Eq. 



V. SPECTRAL DENSITIES AND PROPAGATORS FOR VARIOUS NETWORKS 

Here we indicate four distinct kinds of uncorrelated random networks with qualitatively different asymptotic be- 
haviors of T (0, x) = Tq (x) a,t X ^ oo, where Tq (x) is the solution of Eq. (|2ip at s = 0. At a; = 0, we always have 
To (0) = 1. At a; ^ oo we have Tq (oo) — lima;^oo lims_>o T (s, x) — Q (0). These four types of networks differ from 
each other mainly by a value of the minimum vertex degree. 

1. If the minimum vertex degree > 3, then identically 8 (x) — 0, and Tq [x) exponentially decays to Tq (oo) = 
(see Sec.rOI). 

2. If q„i = 2, then identically 8 {x) = 0, and Tq {x) decays to Tq (oo) = 0, but slower than any exponent (see 
Sec. EH). 
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3. If qm = 1, then there are two possibihties (see Sec. IV Cj) : 



(a) If zi ~ ip'i (1) > 1, then < Q {0} = tc = (fi (tc) < 1, Jo (x) ^ tc a,s x +oo. In this case the graph has a 
giant connected component and a number of finite ones. 

(b) If zi = (p[ (1) < 1, then Q (0) = 1. Tq (x) — 1 a,s x > 0. In this case the graph consists of only finite 
connected components. 

Let us assume qm > 1 and consider in Eq. (j2ip the case of small positive s and large x. According to definition 
P0|) . T {s,x) is actually a Laplace transform of the probability distribution of the non-negative random variable r. 
Hence it cannot decay at a; ^ oo faster than exponentially. In Appendix |E] we show that 



T {s,x) — > A exp [~Tm (s) X ~ (s, x)] , 



(48) 



Here A is simply a constant, and is some correction term in the exponential. The coefficient at the main, linear 
in a;, term in the exponential turns out to be the same as for regular Bethe lattice. It is defined by the relation: 



+ (qm - 1) T„ 



(49) 



This equation has two real solutions as s > Sc, where 



2v^ 



qn 



1 < 0. 



The physical branch of r,„ (s) is the branch, positive at s > 0. The other term in the exponent in Eq. 
sublinear function of x. Namely, 



i9{x) = Bx", 



In {qm - 1) 



In [1/ (1 - Te 



21n[l/(l-r„0] In [1/(1 



0]^ 



(50) 



(51) 



where B is some constant. Here we introduced Tc = Tm (sc) = 1 — ^l\/qm — 1- 

As s is close to Sc, a is close to 1, and t?(s,a;) becomes comparable with the main term. It is this region that 
determines physically interesting results. The behavior of -d (s, a;) at large x and s close to Sc determines the behavior 
of the spectral density p (A) near its edge Ac = — Sc and the behavior of the autocorrelator po (i) at large t. It turns 
out (see Appendix [E|, that the analytic continuation of 'd {s, x) on negative s < —Ac, ^ (—A, x), as a function of x is 
singular at some Xs ~ 1/ (A — Ac). Therefore the upper limit of integration in Eq. (j22p is in the upper half-plane of 
X, Imx > for ImA > and vice versa. Then from Eq. p4)) we obtain 



P(A) 



;P:e^>[T(-A,x)]. 



(52) 



Note that if qm = 2, then Sc = Tc = 0. We consider this case separately in Sec. IV Bl 



A. Minimum degree gm > 2 



Let us set T (s, x) — A exp [—TcX — 1? (s, x)]. Note the difference of the definition of with that in Eq. here 
we have the term —TcX instead of —Tm (s) x in the exponent. Therefore d' (s, x = +00) = Tm (s) — Tc <C 1 now is not 
equal to 0. In Appendix IE II we obtain the following expression, valid when a; ^ 1 and \s\ ^ 1: 



(9m - 1) 

where C ~ 1 is some number. Replace s by —A < —Ac, Ac 



- Sc In (Cy) 



a/4 



In (q„ 



■d (—A, x) = T7T cot 



Sc — 2\/ qm — 1 — 

yX^clnjCy) 

.1/4 



Then we have 



(qm - 1) ln((?m - 1) 



{q,n - if^ 

This function has a singularity when the argument of cot equal to tt, i.e., at x = xq, where 

TT {qm - 1)^^"* In {qm - 1) 



xq = C ^ exp 



VA^ 



{qm - l)-(^'"-l)^'VyA^ 



(53) 



(54) 



(55) 
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When X is close to Xq, one can replace cotz ^ — l/(7r — z) in Eq. ([54)1 . 

Since T (—A, x) is small at large < x < tt, one can replace f [T (—A, x)] by its leading term 11 tp'^'^T (—A, x). 
Then, changing in Eq. ([52|) the integration variable, x = x^y, and taking into account Eq. ([48|) . we obtain up to a 
factor ^ 1: 



, , I dy 

p{\) ^ xo — exp 

' ZTTl 



'Xoy 



a 



qm In {qm - 1) 



1 



6 = 



1 



(56) 



Finally, calculating this integral in the saddle point approximation, we obtain the density of eigenvalues of the 
Laplacian spectrum near its endpoint Ac: 



p (A) exp 



/3 



2VA - Ac 



dexp 



VA^ 



, /? = TT {q,n - 1)^^* In {qm - 1) , 



(57) 



where d is some constant. Substituting Eq. ([57|) into the expression for the autocorrelator ([15)) and using the saddle 
point approximation to calculate the integral, we get: 



Po (t) exp 



-Xrt 



f3H 



W {dt) 



(58) 



Recall the notation Tq{x) — T{Q,x). Since r™ (0) = (q™ — 2) / (g„i — 1) > 0, we have i^'^ [To (x)] ^ 
exp [— ((7m — 2) Tm (0) x] at X ^ +CXD, and the kernel of the integral equation (p9)) satisfies the condition (l30l) . It 



implies that at s = in the discrete sequence of characteristic numbers Pm (0) 



there is the minimum one, 



/iQ > 0. In Appendix IF] we show that (i) /ip = 1, (ii) this characteristic number is the minimum one, and (iii) the 
corresponding normalized eigenfunction is 



■00 (0, x) = -00 (x) = -doTo [x) ) do 



dx ^'i[To(x)]T^2 



-1-1/2 



(59) 



Here do ensures proper normalization (|3ip . and the minus sign stands simply for convenience ensuring ip^ (x) > 0. 

When s > Sc, in particular, near s — > Sc, the kernel in the integral equation ([29]) is well-behaved, and 
all Pm (s) are analytic functions of s. We can leave in Eq. (|34p only the leading term with the minimum p^. Then 
Pi (s) ~ Aq (s) Bq (s) (s) for large distances I from the initial vertex. So at large time t and large distance I, the 
propagator pi (t) is approximately 



Piit) 



+ ZOO 



ds 
27ri 



e^*Pi is) 



-t-200 



ds 
27ri 



e^Mo (s) Bo is) (s) . 



(60) 



If the expression under the integral is analytic in s along the integration contour, the main contribution to the 
asymptotic of the integral gives the vicinity of the saddle point, where st — Hn/iQ (s) is maximal. The saddle point 
position Sc is the solution of the equation t = Ipl^ (s) //iq (s). As a result, we have 



Pi{t) 



1 



x/27r/3 (sc) I 



Ao (sc) Bo (sc) Mo (sc) exp [sc^ - Hn^o (sc)] 



(61) 



where (3 (sc) = [In/^o (•s)]"|^_^ • At a given t 3> 1, this expression has a maximum as a function oil — Im it) at the 
point where 



d_ 
di 



[sct -llnpo (sc)] = In Po (sc) = 0, 



i.e., where po [sc {Im, t)] = 1- Here Sc (Z, t) is defined from the saddle point condition. Since po (0) = 1, the propagator 
pi (t) is maximal at / = /,„: Sc (Irmt) = 0. The behavior of po (s) at small values of |s| determines the shape of the 
propagator near its maximum point. Since po (s) is an analytic function near s = 0, and po (0) 1, one can write 
\npo (s) = as — /3s^/2 + • • • and replace Ao (s) and Bo (s) by Ao (0) and Bo (0). Then the expression ([M]) is reduced 
to a Gaussian integral, and we have 



_ ^0 (0) Bo (0) 



(t-g/)^ 
2/3; 



(62) 
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On the left-hand side of the normahzation condition (|36p , only the term with m = has a simple pole singularity 
at s = 0. Then we have hm^^o (s) Bq (s) / [fiQ (s) — 1] = Aq (0) Bq (0) /a = 1. We substitute the expressions for 
Aq (0) and Bq (0) from Eqs. ([33|l and (|35ll . where the function t/jq (x) is expressed in terms of Tg (x) by using Eq. 
(1551. This leads to 



(O)Bo(O) 



(x) ^[ [To (x)] (1 + ^) r„ J- dx ^' [To jx)] U (^) 
j^dx^'^[To{x)]T;,^{x) 



(63) 



The parameter (3^1 must be positive to ensure the convergence in the summation over I. Equation (j62p . as one 
can sec from its derivation, is valid if the saddle point position |sc| = \{t ~ V")/(/50l we may replace (31 in 

Eq. ((62|) with its value dX I = Im, Pvt, and, finally. 



1 



\/27rDt 



exp 



2Df 



(64) 



where D = (3v^. Despite our network is random, a signal spreads over the network as a Gaussian packet, moving with 

the constant velocity v from an initial vertex, and with the dispersion [l — Im)'^ , which grows linearly with time. This 
is the same kind of evolution as on a regular Bethe lattice. 

Equation (I64p is valid when one can neglect terms of the order of s"^ and higher in the expansion of In/xg (s) in the 
powers of s, i.e., l\s\^^ -^t\s\l < 1. Since Sc = [l/v — t) /I ^ {I — vt) /t, this condition is reduced to \l — vt\ <C t^'^. 
The width of the packet is ~ t^^^ <C t"^^^, and so expression is relevant. 



B. Minimum degree Qm. ~ "2 



If qm = 2, then Sc = as one can see from Eq. (|50p . That is, T{s,x) becomes nonanalytic at s < 0. Besides, 
Tc = T,n (sc) ~ 0, SO that the decay of Tq (x) = T(0, x) is nonexponential in contrast to Qm > 2. Setting T(s,x) = 
exp [—1? (s,x)], we obtain the following expression for small s and large x > (see Appendix IE 2^ : 



i?(s,x) 



1 



\/ sx {a /it + sx) H — arcsinh 

TT 







+ \ln(sA 


--) 




TTX/ 



= vr In 



2n(2) 



>0, 



(65) 



where C 1 is some constant. In the following we omit numerical constants as inessential. When analytically 
continued to s = —A < 0, i9(— A,x) as a function of x acquires a singularity at x = Xc = a/nX. The density of 
Laplacian eigenvalues, p (A) , can be obtained from Eq. ((52|) . The main contribution to the integral in Eq. (|52|) arises 
from the close vicinity of the singularity point. In Eq. (|52p . we expand near Xc in the integral and change the 
integration variable from x to C = A (xc — x). This results in 



P(A) 



1 



A3/2 



exp 



— 

VA/ Jc 27rz 



Here the integral term is ^ A^/^, and the asymptotics at < A <C 1 is 

1 



A4/3 



exp 



Va) 



(66) 



We substitute this expression into Eq. (|15p . and by using the saddle point approximation, arrive at the following long 
t asymptotics for the autocorrelator: 



exp 



(67) 



Let us now consider the propagator pi{t) at I ^ 1, t 3> 1. This asymptotics is also defined by Eq. (pO)) . As for 
Qm > 2, the main contribution to the integral is from the region of small \s\. The difference is that here Aq (s). Bo (s) 
and fio (s) £^11 have a singularity at s = 0. Namely, s = is a branching point, giving a cut along the line (0, oo) in the 
complex plane of the variable s. We will show, however, that this singularity is very weak and does not contribute 
essentially to the propagator, except of relatively small distances 
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Indeed, the small s, large x asymptotics of the eigenfunction i/jq (s,x), corresponding to the largest characteristic 
number {s) = 1 + o (s), is (see Appendix IE 2^ : 



^0 (s, x) w x~^/^T (s, x) - x~^/^ (s + —) 

\ nx / 





ax -1/4 / 






cxp <^ 






TTX/ \ 





\/ sx {a/n + sx) H arcsinh \pKsxfa 



(68) 



Then, comparing the leading terms in Eqs. ([33|) and ([35|) with that in Eq. ([22|) . we conclude that the asymptotics of 
Im Aq (—A) and of ImSo (—A) on A are nearly the same as that of p (A) ~ ImPg (—A). The difference is in powers of 
A in the pre-exponential factors. In the leading order. 



Im[Ao (-A)] ~Im[Bo(-A)] 



^e^-T{-X,x)^o{-X,x) 



A-^/^expf--^ 



(69) 



The rate of singularity of /xq (s), if measured as a jump of a function across the cut near its branching point, is even 
smaller than in Eq. for small A = — s > 0. Let us take the eigenfunction equation ([^5]) at m = 0, setting x — 0. 
Then we have 



Mo ^ (s) V'o (s,0) 



dy e-(i+^)V'i[T(s,y)]Vo (s,y) 



Then, setting s = — A < 0, and properly deforming integration contour, we obtain in the leading order: 



r dx 

ImV'(-A,0) - ^/'(0,0)Ini/xo (-A) - / — exp 

Jc yx 



1 



-(1- A)a; -^{-\x) 

Va 



where the function i? is given by Eq. (|65p . As a function of x this integral has a singularity at a; = o/ttA. In comparison 
with the integral for p (A), the above integral has an additional term —x in the exponent, which turns into —a/-K\ at 
the singularity point. Therefore, we estimate the singularity of pq near s = as Im/ig (^A) ~ exp (— o/ttA). 

Since all multipliers in Eq. (j60p have sufficiently weak singularities, we replace Aq (s) and Bq (s) with their values at 
s = and neglect the singular part of /iq (s)i leaving only the regular part of the expansion: In/ip (s) = s/v+Ps'^ /2+- ■ ■ . 
As a result, we arrive at the same Gaussian expression for the propagator, Eq. (j64p . 

If we, however, fix the distance I ^ 1 and increase the time t, the saddle point Sc < in the integral ([60]) moves 
farther in the direction of negative s, and at large enough t the contribution of the singularity becomes essential. 
Deforming contour of integration, we rewrite Eq. (|60p in the following form: 



Piit) 



dX 



-At 



Im [Aoi-X)Bo i-X)pl-' (-A)] . 



(70) 



Im/iQ ~ exp (— o/ttA) is small compared to Im Aq (^A) ~ Im_Bo (^A) 
of pq and set Inpo (—A) = —X/v. Thus we arrive at 



Piit) 



dX 



exp 



-Alt-- 

V 



exp(— a/\/A). So we neglect the singularity 



Calculating the integral in the saddle point approximation we obtain 



Pi it) 



I 

t 

V 



-5/18 



exp 



-3(^ 



2/3 



I 

t 

V 



1/3 



(72) 



Expanding hipo (s), we neglected terms of the order of and higher. This is justified if the saddle point position 

in the integral ((7T|) . A^ ~ (t — l/v)~'^^^, obeys the condition IX^ ^ Xs ^^"^ which is equivalent to t — Ijv ^ t^l'^. 
Otherwise, pi (t) is given by Eq. (j64p . which means that the probability for the signal to return is small. This form 
of the packet tail is due to the possibility that either initial vertex or the final one in the l-th shell may occur in a 
chain fragment in the graph. 



C. Minimum degree = 1 

When there is a finite fraction of "dead ends", i.e., vertices of degree 1, the network contains finite-size connected 
components. They lead to the (5-functional peak in the Laplace spectrum and so to nonzero limits of the averaged 
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FIG. 4: Graphical solution of the equations for tc, tc ~ f i (tc) and for ts > tc- (fi'i (ts) = ipi its) /tg. 



propagators at t ^ oo. If (p[ (1) > 1 f Appendix [B]) . then besides the finite connected components, there is a giant 
connected one whose size scales as the network size. Here we show that the contribution of this giant connected 
component to the observable quantities is qualitatively the same as in networks with = 2. 

If q-m = 1, Eq- (EB) still has the nontrivial solution Tq (x) = T{0,x). To (0) = 1 as for any other Qm, but 
To = lima;^o hms^o T {s, x/s) = 9 (0) =tc>Q (see Appendix |B]). At small s > and large x > 0, the function 

T (s, x) is close to Q (sx), and so we search for T (s, x) in the following form: 

T{s,x) = e{sx)+e-^^'-''\ (73) 

where the last term is assumed to be small. The asymptotic solution for i) is (Appendix IE 3[) 



^{s,x) = ^g{sx) + ]\nsg'\sx)+C, g {z) = f dyJl -^^^El^M ^ (74) 

where C ~ 1 is some constant. Continuing this result to s = —A > 0, we take into account that g (2) has a singularity 
at z = < 0, where Zg satisfies the equation 1 — \a.ip'i[Q{zs)]/ Zg — 0. The equation for Zg, ip'^ [8 (2^)] — becomes 
more comprehensive with a new variable tg = Q (e~^=). Using the implicit definition ()38|) of O (2), we arrive at the 
equation for tg: 

[ts) = — . (75) 

This equation is shown graphically in Fig. |31 together with the equation for tc = fi (tc) < tg, Q (tc) = 0. 

An expression for the spectral density p (A) may be obtained by calculating the integral in Eq. ([5^ , where for small 
A and large x, we approximately set: 

if [T (-A, x)] « ^ [6 (-Ax)] + [9 (-Ax)] g-''^^^'"). (76) 

While the first term results in the (5-functional peak, the integration of the last one gives the asymptotics of p (A) 
at small positive A. The main contribution into the integral gives the close vicinity of the positive singularity point 
Xg = —Zg/X > 0. Near this point the first multiplier in the last term in Eq. (j76l) can be replaced by a constant 
if' [Q (2s)] = (p' (tg), and the function g in the expression for z9 (|74|) can be replaced by its expansion g (zg + t]) ^ 
ia + ihif/'^, 1, |?7| < 1 with 



*-^^ ^i(x)-V,(x) / y. i-) . (77) 
xipi (x) y ln[(/?i (x) /xj 

(see Appendix IE 3|1 . Then we arrive at the following asymptotic result for p(A): 

(78) 



p(A)-p(^'^)^(A)^A-^/i"exp {-^ 
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As it follows from Eqs. ((78)) and ([15]), the autocorrelator pq {t) decays to its equilibrium value as 



(cq) 
Pq 



dX 



exp 



-Xt 



t 



-7/30 



exp 



2/3 



1/3 



(79) 



One can calculate the propagator (t) at large I using Eq. ([50)) . As compared with Sees. IV A) and FVB) the kernel 
of Eq. ([^5)) is not any more bounded at s = because the integral in Eq. ([50)) becomes divergent. Due to this fact, the 
spectrum of Eq. ([29]) contains continuous part. Let us find eigenvalues /zj and eigenfunctions ■0^ (x) in the continuous 
spectrum. (The notation V'm {x) = tpm {s — 0, x) we leave for the discrete part of the spectrum.) We saw in Sec. IIVI 
that when s ^ 0, the recursion relation ([^ can be transformed to Eq. ([i5)) . assuming that Ui {s,x) « ui (sx) at 
small s. An equation for the eigenfunctions is 

i^i (x) = Mee-'^t^'i [6 (x)] (x) , 

which has the solutions ipi^ {x) = (5 (x — ^) corresponding to the eigenvalues — /f'l [6 {x)]. It is the continuous 
part of the spectrum that after proper modification of the relations pm35p , gives the stationary part of propagator 
()46p . Suppose that there is a giant connected component in the network. Then along with the continuous part of 
the spectrum, whose minimum characteristic number is ^J-^=o = ^/fi (^c) > 1 = /^m=05 there is a discrete spectrum 
with the minimum characteristic number /iq (s = 0) = 1 corresponding to the eigenfunction (x) = — Iq (x) (see 
Appendix [F| . 

In the same way as for q„i ~ 2 (see Appendix lE 3p . one can show that the asymptotics at small s and large positive x 



of the eigenfunction i/jq (s, x), corresponding to the lowest eigenvalue fxo (s), is tpo (s, x) ^ x 
'd (s, x) is given by Eq. ([7i)) . 

From Eqs. ([55]) and ([55]) we obtain 



-1/2 



exp [— ?? (s, a;)], where 



Im [^0 (-A)] ^ Im [Bo (-A)] ^ 

J-ioo 27^^ 



re^^^oi-X,x) 



-|-ioo 



dx 



X 



-1/2 



exp 



(80) 



This equation differs from Eq. (I69p . because in the singularity point Xs = —Zg/X ^ 1, T(— A,Xs) ~ O (— ^s) ^ 1- So 
we omitted T {—X^ Xs) in Eq. (|80p . in contrast to the case = 2, where the function T(— A,Xs) has is of the same 
order of smallness as ipo (— A,Xs). Here, as for qm — 2, the singularity of /ig (s) is such that the jump along the cut 
(— cx), 0) in the complex planes s behaves as Im fiQ (—A) ^ exp (—a/A). 

The derivation of pi (t) for = 1 is similar to that for = 2. We arrive at the same moving Gaussian packet (|64p . 
The only difference is that now we must take into account the contribution of finite clusters (continuous spectrum). 
The results for I ^ 1 and t ^ 1 are 



1 



V2TrDt 



exp 



ji - vty 

2Dt 



(81) 



for \vt — Z| ^ t^/^, where v is given by Eq. ([T7)). and D — v^fi^ (0). In the low / tail, I < vt, vt — I <^/^, the form 
of the propagator is modified to 



PI it) -Pi 



(cq) 



I 

t 

V 



-13/30 



exp 



1/3 



(82) 



Thus, again, we have the Gaussian packet, Eq. 
is supplied with a small tail at 1 <C / ^ i, Eq. 
"dead ends" . 



moving within the giant connected component. This Gaussian 
The reason for this tail is that initial or final vertices may be 



VI. SUMMARY, DISCUSSION, AND CONCLUSIONS 

In this article we have presented a theory, which enables us the analytical calculation of statistical properties of the 
Laplacian operators of infinite random networks and random walks on them. We have considered the resolvent of the 
Laplacian and the propagator of a random walk. These characteristics are connected through a Laplace transform, 
Eqs. ()12p and ([TS]) . In particular, the average values of the diagonal element of the resolvent matrix give us the 
spectral density of the Laplacian, Eq. ([H]) . and the time dependence of the autocorrelator. We have also derived 
equations, which solution allows us to find the averages of the nondiagonal elements of the resolvent. After the Laplace 
transformation, these averages show how the distance of the signal from its origin changes with time, Eq. (Ilip . 
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Our scheme is based on equations relating the distributions (or other statistical properties) of random variables. 
This is an essential advantage over most of existing approaches, based on equations relating the values of some random 
variables for a given network realization. To solve the problems of the Laplacian spectrum and of random hopping 
motion, one must make the following steps. 

(i) Solve the integral equation pT|) for the function T {s, x) defined by Eqs. (|19p and ([20]) . [In the equivalent form, 
it is Eq. (jC7|) .] Technically, it is the most difficult step. We have only obtained the asymptotics oiT{s,x) at 
Rex — > +00. We have found that T as a function of x is an analytic and exponentially decaying function as 
s > Sc, where Sc = — Ac < is a parameter which depends only on the minimum vertex degree Qm- 

(ii) With T{s,x), one can (a) calculate the average of the resolvent's diagonal, Eq. (j^ . then (b) analytically 
continue the result from the positive s to s = — A±iO, A > 0, and finally (c) obtain, using Eq. HH), the spectral 
density of the Laplacian [2^ . 

(iii) With the known p{X) near the spectrum edge, obtain the asymptotics of the autocorrelator pQ (t) at i — > 00 by 
using Eq. 

(iv) Find the sequence of functions Ui {s,x), I > 1, definition ([24| . by using the integral recursive relation (f25|) with 
the initial condition (|26p . Then obtain the Laplace-transformed propagator Pi (s) by calculating the integral 

(v) Calculate the inverse Laplace transform of Pi (s), that is, the propagator pi (t). The asymptotics of pi {t) at 
large I and large t is determined by the smallest characteristic number fj,o (s) at small |s|. 

The results of these calculations of asymptotics are summarized in Table |T] and Fig. [1] If > 3, the tail in the 
density of eigenvalues decreases extremely rapidly with 1/ (A — Ac), see Eq. ([571) . and therefore practically cannot be 
revealed by numerical methods. Studies based on these methods usually result in a form of p{X) resembling Wigner's 
semi-circle law (see, e.g., Refs. [1, §)• This is also the case in networks with q„i ~ 1, 2. 

When are our analytical results observable? Let us inspect the resulting expressions for the propagator pi (t) . Our 
results are based on the tree ansatz: pi (t) should have nonzero values in the small (compared to the whole network) 
vicinity of the starting vertex 0, so that we can treat this region as a tree. At large times t the signal spreads at the 
distance I = vt ^ t, Eq. (|64p . The mean intervertex distance in the network is In [2^ [2^. So, our results are 
applicable if 1 ^ t ^ InA^. In networks with > 3 the decay of the autocorrelator is basically exponential with 
some correction [see Eq. ([55]) ]. This correction can be observed if 

1/ln^i- l/ln^lnA^^ 1. (83) 

It seems to be impossible to fulfil this criterion either in real-world networks or in numerical simulations. 

In the networks, containing chain- like segments, i.e., when = 1 or 2, the criterion is much less stringent. We 
require that the value of the autocorrelator po (t) (|57)) at the characteristic time t ~ InA^, essentially exceed its 
equilibrium value pQ{t = 00) ^ 1/N for a finite network. So in these networks, our dependences are observable if 

ti/3/lniV-ln"2/^iV< 1, (84) 

which is much easier to satisfy than condition (|83p . 

In many applications of the Laplacian spectrum, results, obtained in the infinite network limit, are of little use. A 
good example is synchronization j8l.[28j. In this problem the lowest, size-dependent eigenvalue of the Laplacian plays 
a key role. Let us briefly discuss the role of this eigenvalue in application to our problems. The process of a signal 
spread over the network consists of two distinct stages. We discussed the first one. In the second stage, the essence of 
the process is the relaxation to the homogeneous distribution, where the probability to find a signal at any vertex is 
the same, namely, 1/A^. In this last stage, \pio (t) — 1/N\ <^ 1/N. In this situation loops must be taken into account. 
Furthermore, in this stage, the knowledge of the Laplacian spectral density is not sufficient. Rather, one should ask: 
what is the probability distribution of A2 (the lowest nonzero eigenvalue)? 

We showed that in infinite networks with minimum vertex degree qm > 2, the density p (A) — for < A < Ac. In 
contrast, in finite networks, Laplacian eigenvalues A^ exist in this range, though only a very small fraction of the total 
number of the eigenvalues. The statistics of this part of the spectrum determines the second stage of the evolution of 
Pio it) to the equilibrium. We believe that this stage may be described in the framework of an approach developed in 
Ref. '23l| for calculation of intervertex distance distributions. We leave this challenging problem for future study. 

In summary, we have strictly shown that the region of low eigenvalues in the Laplacian spectra of uncorrelated 
complex networks and the asymptotics of random walks on them are essentially determined by the lowest vertex 
degree in a network. 
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APPENDIX A: OTHER LAPLACIANS AND PROCESSES 



Three different forms of a Laplacian operator are discussed in literature. In this paper we discussed the form ([T]) 
corresponding to the process defined by Eq. ([2|). The second form, 

Lij — Sij ^ijj (-^1) 

corresponds to the following process: 

^ 1 

Pij (0 = —AikPk]{t) - Pij{t), Pij{Q) = Sij. (A2) 

k—l 

This is a random walk process with the unit escape rate of a particle from any vertex. The particle jumps to any of 
Qi nearest neighbors of vertex i with the same probability 1/ft. We do not consider this process here, although it can 
be described in the framework of the approach of this article. We have found that the singularity of the spectrum at 
the lowest eigenvalue of this Laplacian and the long-time asymptotics of the autocorrelator of this random walk are 
quite similar to those we found for the operator ^ and the process ©. 
The third, "normalized" , form, 

Liij — 6ij Aij (-^3) 

(see, e.g., Ref. Q) is, one may say, equivalent to the form (jA2[) in the following sense. Operators (|A2p and (|A3p are 
connected by a similarity transformation. The connecting operator W is diagonal: Wij = SijQi. These two operators 
have the same spectrum of eigenvalues. Their eigenf unctions are connected by the operator W. 



APPENDIX B: DEGREE DISTRIBUTION IN Z-REPRESENTATION 



The Z-representation of a discrete random variable qi = 0, 1, 2, . . . is defined as 

Af oo 



(Bl) 



9=0 



Lp{z) is also called the generating function of n(q). It is obvious that (p(l) = 1. Differentiating (p{z) and setting 
z = 1, we obtain an expression for the average vertex degree, 



^'(l) = ^gn(5)=.(g)^q-. 

q=0 



In general. 



dz 



ip{z) 



For branching numbers bi — Qi ~ I we have 

N 

J- 



AT 1 -'^ -, oo 

1 = i E (^^.-'0 = ^ E ('^.-"-^) = J E (q) 



9=0 



(B2) 



(B3) 



(B4) 
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The function Lpi also obeys a normalization condition, i^i (1) = 1. 

This function was successfully used by Newman, Strogatz and Watts [1^ (compare with the earlier works by MoUoy 
and Reed, Ref. (s^jllil) in their calculations of the size distributions of n-th connected components of a vertex. Recall 
that this is a number of vertices which are not further than n steps from a vertex. For example, the distribution for 
the first connected component in Z-representation is zip (z), for the second one, it is zp [zipi (z)], and, in general, the 
distribution for an n-th component is Gn {z) — zip [Hn (z)]. Here the sequence Hn is defined by the recursion relation 

i7„ (z) = zpi [i7„_i (z)] , Hq (z) = z. (B5) 

Its stationary solution H (z) satisfies the equation: 

H (z) = zvji [H (z)] . (B6) 

So G (z) = z(/9 [H (z)] = l^z^'^^) is the transformed probability function that a randomly chosen vertex is in a connected 
component of size Mi. 

The function G(z) allows one to find, in particular, the relative size of a giant connected component, moo — Nao/N . 
Let us consider the solutions of Eq. (jB6p as z ^ +0, H (+0) = tc] tc ~ (pi (tc). Beside the trivial solution equal to 
zero, there is another solution, tc < I (see Fig. S]): 

tc = m (tc) , < < 1, ifzi= f'l (1) = > 1. (B7) 

ip' (1) 

So ip (tc) is the total relative size of all connected components of the network, and the relative size of the giant 
connected component is 

m^^l~ip {tc) . (B8) 

Note that the condition (jB7p may be written as 

^g(g-2)n(g) >0. 

If there are no "dead ends" in the network, then tc — p> {tc) — 0, and almost all vertices in the network are in the 
giant connected component. 

APPENDIX C: EQUATION FOR THE DISTRIBUTION OF r AND AUTOCORRELATOR 

If n > 1, Eq. ITBl) may be written as (see Fig. [3]) 

sP„+i,, (s) - [P„,j (s) - P„+i,, (s)] + J2 [^"+1.^ («) - ^"+2,fc (s)] = 0. (CI) 
Dividing both parts of the equation by Pn+iA (s), and taking into account the definition (|19p . we get 

'''n,ij (s) 



1 Tn,ij (s) 

If n = 0, Eq. takes the form: 



sPo (s) + J2 - ^ifc («)] = 1- (C3) 



k=l 



Dividing both sides of Eq. (jC3|) by Pq, and taking into account Eq. (fT9|) . we obtain 



90 



S + ^ To,ofc 



(C4) 
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Recursive relations (IC2p express the set of random variables Tn,ij, n > 0, in terms of the set of independent and 
statistically equivalent random variables qm.iy m > n. It is important that the variable t^jj depends only on the 
degrees of vertices belonging to the tree branch, which grows from the edge (n, i) — {n + So in Eq. ([6]), Tn^ij (s) 
is expressed through qn+i,j independent random variables: the branching number bn+ij = Qn+i,j — 1 and &n+ij 
statistically equivalent variables Tn+ijk, k = 1, . . .bn+ij- In the thermodynamic limit, the statistical properties of 
branches, starting at any distance from the initial vertex, are the same. Consequently, all random variables Tn.ij are 
distributed equally, independently of i,j and n. Then, omitting unnecessary indices, one can rewrite Eq. (IC2|) as 



exp 



y 



l-r(s) 



b 



Yle-y'-'^'l (C5) 



-sy 

k=l 



The next step is averaging both the parts of Eq. (|C4ll . We use definition ((20|l , properties of statistical equivalence of 
T, and mutual independence of the branching number b and all r^,. We also use the following integral identity: 

-y/c Vy dx 



g-y/o ^ Vil , r^^ (2V^) e"^ (C6) 



ioc+S V 



where Ki is the MacDonald function of index 1. Then 

y \ ^ VI 



Finally, we have 



^yVv ^K,{2^)e^T{s,x)^e-^y({[T{s,y)]'})^e-^y^^[T{s,y)], (C7) 

where definition (jB4p was used. (Here 6 is a branching coefficient of some edge.) 

Introducing ^ (s) = r (s) / [1 — t(s)], so that r (s) = (s) / [1 + ^ (s)], enables us to use the following integral 
identity equivalent to Eq. (jC6P : 

e-/" = 1 + ^ r ^I, (2^) e-^y. (C8) 

Jo \/y 

Here Ii is the modified Bessel function of index 1. So 



'<^'^) = l + V^J^ ^h{2^)> 
Note Eq. (|C2p leads to the relation: 

6 

as)^s + Y.^k{s). (C9) 



k=l 



Averaging in the integral in the same way as in Eq. (jC7p . we arrive at Eq. (j2ip . 
After the averaging, Eq. (|C4p takes the form: 



<.oo / 90 

Po (s) = / dx e""^ ( TT exp [-xtq.o/c] 

•^0 \fc=i 



Taking into account the property of statistical independence and equivalence indicated above, we obtain Eq. 



APPENDIX D: DERIVATION OF THE RECURSION RELATION 



Let us consider the following expression: 



is) 



exp 



L '^n,ij (^) 



(Dl) 
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Differentiating the identity (|C6p with respect to x we have 

-e"^/"-/ — ifo(2V^)«e"^ (D2) 

Substituting a = 1/ (1 — r) and using the definition of [/;, Eq. ((24)) . we transform the expression IjDip into 

eM — Xo(2V^)5„,„exp{x[l-v,,(s)]}, (D3) 

J~ioc+5 

where ii'o is MacDonald's function of index 0. In the infinite network the expression (jD3l) is independent of the chosen 
edge (n^i) — (n + and depends only on / and on s. Averaging Eq. (|D3[) we get 

+ioc+S 1 

—Ko{2^)e^Ui{s,x). (D4) 

-ioo+5 



On the other hand, due to the tree-hke structure, the first multipher in the angular brackets in Eq. (jDl|) may be 
expressed as a sum of terms with I ^ I — 1 (see Fig. [3]): 

i''^) = [1 - ^n,, is)] (s) . (D5) 

fe=l 

Using Eq. ([5]) together with Eq. (|C2p . we see that the expression (|Dip is equal to the following one: 

fc— 1 m— 1 

Let us average Eq. (ID6|) taking into account the statistical properties of the variables r and b (or g), indicated above. 
Note that in each of ^n+ij terms we have bn+ij — 1 multipliers (exp[— yr„+i (s)]) — T(s,y) with fc 7^ m, and 

the multiplier (<S'^'_(_i (s) exp[— yr„+i (s)]) = C/j-i (s,y). So the remaining average over fon+ij = ^ can be easily 



■J 

performed, which gives 

e-^y J2 b [T {s, y)]'-' ni(6)[7z.i(s, y) = g-^V'i [T (s, y)] Vi-i (s, y) . (D7) 

6 

[Recall that the distribution function of h is Hi (6), Eq. (fT7|) . i.e., (z) in Z-representation, Eq. (jB4p .] Equating 
expression (jD4p to Eq. (fTT|) . that is a different representation of expression (jDl|l . we derive the recursion relation for 
[/„ in the following form: 

r.+ioo+i5 J 

— (2V5y) e-C/, (s, x) = e-^V; [T (s, y)] C/,_i (s, y) . (D8) 

-ioo+i5 

Equation for eigenfunctions il)ra (s) can also be written as 

— (2Vi^) (s, 2:) = e-^ Vi [T (s, y)] (s, y) . (D9) 

-ioo+(5 ""^ 

Let us now replace t in the definition of [//, Eq. (|24p . with its expression in terms of the random variable ^ (s), 
r (s) = ^ (s) / [1 + C (s)]. In turn, for ^ (s) we use relation (|C9p . Again, use Eq. (jPSp for S*. Differentiating integral 
identity (|C8|) with respect to a; gives 



/•OO 

/ dyIo{2y^)ae-"y. (DIG) 
Jo 



Then we have 



•^0 \ fc=i ™=i 

= e-^/ dy/o(2Viy)e-(i+^)^ V6[r(s,y)]^ni (&)[/,_! (s,y). (DU) 
Jo 

Using Eq. (jD7p for averaging over b readily leads to Eq. The initial condition (1^^ follows directly from the 

definition of Ui, Eq. 



21 



APPENDIX E: ASYMPTOTIC SOLUTIONS OF INTEGRAL EQUATIONS 

Calculating the asymptotics at large x we replace the MacDonald functions [2,Jxy) with the leading term of 
its asymptotic expression: 



(El) 



{z) - 

This asymptotics is independent of v. Then Eqs. (|C7p and (jD8|) at large x take the forms: 

■ exp [x - 2^/xy) T [s, x) = e'^^^ipi [T (s, y)] 



yjj^ /•+'°°+^ dx 

3/4 



and 



1 



100 + 5 



Equation (|D9p in the asymptotic limit has the form: 
1 



zv27r J-ioc+s ^ 



+ioo+S 7 

— exp {x - 2^) Ui {s, x) = e-^Vi [T (s, y)] Ui^i {s, y) . 



+joo+<5 

Yyj exp (x - 27iy) V'm (s, x) = (s) e'^^v^i [T [s, y)] iprn (s, y) 



(E2) 



(E3) 



(E4) 



Let us first consider Eq. (jE2| in the case s > 0, y — > oo. According to the definition ([20|) of T {s,x), this function 
is the Laplace transform of the probability density of a random variable t (s). This variable satisfies the condition 
< r < 1, as it follows e.g., from the recursion relation (jC2|) . Hence, (i) the function T{s,x) is analytic everywhere 
in the complex plane x, and T (s, x) ^ as Re a; — > oo; (ii) T (s, x) cannot decrease with x faster than exponentially. 
Then T (s, x) can represented as 



T{s,x) ^ A exp [-Tm (s) X - z9o (s, x)] 



(E5) 



where r„i > 0, and dx'^^o is,x) — > as a; ^ +oo. If q„i > 1 (we will consider this case separately), then ipi on the 
right-hand side of Eq. (jE2p can be replaced with its leading term, (pi (z) {<lm) z^"^^^: 

9~172 / (^^P[{^-Tm)x - 2^/xy - 'doix)] = -^n(g,„)yl'''" exp[-~{l+s)y - {qni-'^)TmV - (gm-l)i?o(2;)] • 

The integral on the left-hand side may be treated in the saddle point approximation. The saddle point equation is 
the condition that the derivative of the function in the exponent becomes equal zero, namely, 

y = x[l^ r„, - {x)f . (E7) 
This equation also expresses y in terms of x. On the right-hand side, we assume that 

^^0 (y) ~ ^0 [(1 - Trnf x] - 2 (1 - Tm) X^'^ (x) [(1 - T^f X 

One must prove afterwards that the neglected terms of the order of xi?[f and with higher derivatives are small. If 
we set = in the pre-exponential factor of the saddle point approximation, it reduces to 1. So we arrive at the 
following equation for t^o^ 



In 



-Il{qm)A''" 



+ (l-r™)2 



S + (gm-l)T„ 



1 - T„ 



X - 2(l-r„) 



S + ((7m-l)T„ 



1-T„ 



x^^(x) 



+ (q™-l)i9o[(l - T„)2x] - Mx) + [l + s + (q„-l)r™]x<(x) - 2(q™-l)(l - T^)xd'„{xy„[{l - t„)2x] = 0. (E8) 
The main term of this equation, linear in x, reduces to zero if 



1 - T„ 



s + {q,n - 1)t™. 



(E9) 
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It also reduces the third term in Eq. (|E8p to zero. Suppose that our network is a the regular Bethe lattice with the 
coordination number Then Eq. (jC7[) . equivalent to Eq. (|22p . has the exact solution T(s,x) — exp [— (s) a;], 
where r™ (s) is the proper solution of Eq. (|E9p . This r™ (s) is a regular function of s as s > Sc = — 9m + "^Vlm — 1, 
Sc < 0, and (sc) = Tc — 1 — l/V'Zm — 1- At s = Sc, Tm (s) has a square root singularity. So in the regular Bethe 
lattice, the density of Laplacian eigenvalues p (A) is nonzero at A > Ac = — Sc = Qm — 2\/9m — 1, and p (A) ~ V-^ ~ -^c 
at A — Ac <C 1. Thus, we can conclude, that for any network with qm > 1, the edge of the spectrum is Ac > 0. 
Moreover, Xc{qm > 2) > 0. In random networks, the asymptotics of p(X) turn out to be sharply different from a 
regular Bethe lattice. 

Requiring that the main correction to the leading term in Eq. (|E9p also asymptotically vanish gives 
((7„j — 1) iJo (1 ^ Tjn)'^ X ~ "do (x). This equality is satisfied when 

, f . R " In(gm-l) ln[l/(l-Tc)] 

Wo (x) = Bx , a — — — p— "7 -r = - — - — -— ; — . (ElO) 

' ' 21n[l/(l-T„0] ln[l/(l-r„)] ^ ' 

If s > Sc, then > Tc and a < 1. This means that all approximations made during the derivation are justified. 
Therefore the integral in Eq. (|22[) is convergent, and Pq (s) is a regular function of s. If, however, s — > Sc, then a — > 1. 
That is, the last two terms in Eq. (jESP should be also taken into account when s is close to Sc- In this region of s, 
Eq. (|E8|) and a similar equation for the asymptotics of tpQ [x) , which can be derived from Eq. (jD9|) , must be treated 
in different ways for > 2 and for q„i — 2. Note that if = 1, then Eq. (|E8[) must be replaced with a slightly 
different equation. 



1. Minimum degree > 2 

In this case we can set the first term in Eq. (jE8[) to 0, properly choosing the value of the constant A in Eq. (|E5[) . 
We set T (s, a;) = exp [— ((5, a;)]. Here i9 ((5, x) — ■dti + t^x includes, besides t^o, also a slowly varying linear term t^x. 
Here we introduce a small variable 5 — ^ s — s^. 

First, let us consider Eq. (jE8p at s = Sc and t = Tc. We have 



= 0. 



V 9m -L / \ 9m ^ / 



Now we make the substitution: d (x) — {q^ — 1) a;x(liia;). We assume that x is a small and slowly varying 
function of its argument. Then we make the following approximations, which must be justified afterwards. Replace 
X [z — In (9,„ — 1)] in the first term with x (z) — x' (z) In {q,n — 1), where z = Iny, and neglect all derivatives of x in 
the last two terms. As a result we get 

ln(9„,-l)x' (z) + x'(^) = 0. 

This equation has the solution: x (z) = In (9,„ — 1) / (z + c), where c ^ 1 is some constant of integration. Thus, 
finally, we obtain 



x In {q„ 



I) 



^/q,n - 1 In (Cx) 



(Ell) 



Now assume \S\^ = \s — Sc = \X — Xc\\ ^ 1- The first term in Eq. (IE8p reduces to S'^y/ (9™ — 1). We neglect the 
second term of the equation, assuming it to be small. After the same set of substitutions and approximations as in 
the case s = Sc, we have the following equation for x (z): 

\n{q,n-l)x'iz)+xHz) = 



Solving this equation, we obtain the following result for d {S, x) : 



((5, x) 



Sx 



coth 



V9^ 



^ln(Ca;) 



(9m - 1)^^^ In (9m - 1) 



(E12) 



After substitution 5 — ^ s — Sc, this turns into Eq. 



23 



2. Minimum degree qm ~ 2 

Here (s) — > as s — > 0. Then at small \s\ we can consider (s, x) = i?o (s, 2:) + '''m (s) a; as a slowly varying 
function. When calculating the integral in Eq. (|E2p in the saddle point approximation, we take also into account the 
pre-exponential factor as a correction, though it is close to 1. Replacing on the right-hand side ipi (T) with its leading 
term, linear on T, and taking into account the saddle point equation y — x[l — 1!}' (x)]^, we have 



1 - (x) 



1-t}' {x) - 2xd" [x) 



1/2 



exp [x] - ^ {x)] 



2n (2) 



exp [-SX - -i? (x) + 2x1?'^ (a;)] 



(E13) 



Here we omitted negligibly small terms: 2sx'd' {x) and others. 

Accounting for the smallness of "d' and xd"^ we obtain the equation: 



xd'^ {x) - yd" {x) = sx 



TT In 



2H (2) 



> 0. 



(E14) 



We assume that the first term on the left-hand side is small and search for the solution of this equation in the form 
■d = -di -\- ^?2- Here di must be found from xd'i = sx + o/tt. At s = we find di — ^^axj-n -I- c, where C is some 
constant of integration, C ^ 1. At s ^ 0, performing the integration, we have 



1 



i9i (x) — —=] {sx) + C, / (z) = J z (a/iT + z) H — arcsinh W — . 

J S TT V Q 



(E15) 



In principle, here C — C (s), but for small s one can set C [s) — C — C (0). For §2 we have 2i9'ji?2 = Therefore 
up to the constant, i?2 = (Ini^i) /2. As a result, we have asymptotically the expression for §^-01+ §2- 

We replace in Eq. (ID9p Kq with its asymptotic (|Eip at large values of argument. Then, taking into account that 
T (s, x) is small at large x, we replace Lpi (T) on the right-hand side with its value of zero argument, 2n (2) /q. As a 
result, we arrive at the following equation for ipo {s, x) = exp [— >f (s, x)]: 



2iy^/\ 



exp [-^^/xy + x~ >t {x)] = fio (s) -H (2) exp [-sy - k (y)] . 

-ioc+5 2; ' q 



(E16) 



This equation differs from Eq. (jE6|) for ?9 (x) = i?o (x) + Tm (s) x only in the pre-exponential factor on the left-hand 
side. Quite analogously to Eq. (|E14p . we obtain an equation for x, 



xyi"^ (x) — xk" (x) — K [x) = sy -\ , a' = tt In 

TT 



Mo (g)g 
2H(2) 



a > 0. 



(E17) 



Here b is given in Eq. (jE14p . Comparing the above equation with Eq. (IE14p . we conclude that x (x) — i?i (x) + X2 (x), 
where i9i (x) is given by Eq. (jElSp . and X2 must be found from 2x{}'i (x) x'2 (^) — ^^1 iv) + ''^i i^)- The solution is 
X2 = {lnx'd[) /2. As a result, accounting for Eq. (jElSp . we obtain the expression (ISS)) for = exp {—x). 



3. Minimum degree qm ~ 1 

To obtain equation for Q (z) = Imis^o T {s , z / s) , let us start with Eq. (|E2p . which is valid as y — > 00. Let us 
replace y in this equation with z/s simultaneously changing the integration variable: x = Q/s. Then we have 




ioo+S 



■ exp 



(x/i-x/C)' 



r( s,^ ) =e->i [T{s,z/s)] 



(E18) 



In the limit s — > the saddle point approximation becomes exact, with the saddle point condition simply = z. 
So, assuming that the limit ([57)1 of the function T exists, we immediately arrive at Eq. ([55)1 for 8 {z). The recursion 
relation (j43p is obtained in the same way by using the asymptotic equation (jESp . 
We can reasonably assume that at small s and large x, 



T (s, x) — Q (sx) + exp [—d (s, x)] , 



(E19) 
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where the last term is small. Substituting this into Eq. (|E2p and linearizing the right-hand side with respect to e ^ , 
we obtain 

/ exp [~2^ + x^d {x)] = v^'i [0 (sy)] exp [sy ~ d {y)] . (E20) 

Repeating the steps leading to Eq. (IE14|) . we get 

xd'"^ {x) - xd" {x) ^ sx- In ip'^ [9 [sx)] . (E21) 

Treating the second term on the left-hand side as a perturbation, we set ?? = i?! -I- i?2, |^^2| ^ The equations for 
di and i?2 are 

xd'^ [x] = sx - In (^'i [6 [sx]] , 2'd[ (x) (x) = i^'/ (x) . 

Then one can easily obtain 



(s, x) = -^g (sx) + C, g (z) = dC^ 1 - ^^2 (s, x) = - In [sg'' {sx)\ , (E22) 

where C ~ 1. Finally, for (x) — -di (x) + i?2 (x) we get formula ([7^ . 

The function (; (z) has a singularity point at z — Zg < 0, when the expression in the square root under the integral 
becomes 0, i.e. when ip[ [9 (z^)] = e^' . Here Zg is defined as = 6 (e~^=), tg is the solution of Eq. ([75]) [see Fig. 

It is obvious that g (zg) = ia is imaginary, because the expression in the square root in Eq. (jE22p is negative. 
The calculation of a may be simplified if we replace the integration variable C with ^ = Q {(,). We use the definition 
of the function 9, Eq. from which C = In VPijX) H] follows, so dC, = [ip'-^ (C) /<^i (C) - 1/^] d^. We substitute 

these relations into the integral for g {zg) in Eq. (jE22|, and take into account that ^ = tc = 9 (0) at the lower limit 
of integration, C = while = tg on the upper limit. This gives Eq. ([77]) . Note that z — tg \& a. singularity point of 
the function 9 (z). Since the derivative of the reverse function z = In [(fi (9) /9] is zero, Q {zg + rj) = tg + O {^/v) 
small \r]\. Therefore in Eq. (IE22P In ip[ [9 {zg -f 77)] / (z^ + 77) - 1 ^ and so g {zg + r]) ^ ia + O (77^/^) . 

The calculation of the asymptotics of the eigenfunction ipQ (s, x) is quite similar to that for = 2. Let us represent 
ipQ (s, x) = exp [—XT (s, x)]. From Eq. (jD9|) we obtain the integral equation for x which differs from Eq. (jE16[) only in 
that the constant (p[ (0) = 211 (2) /q should be replaced with the function (p[ [9 (sy)]. Proceeding further, we have 

xk''^ (x) — xh" (x) — x' (x) = sy — In ip'-^ [9 (sx)] . 

This equation for x differs from Eq. (IEI7P only by the last term on the right-hand side. So, as at q^ — 2, we have 
asymptotically ■00 {s, x) ^ x^^^^T (s, x) ^ x^^l'^ exp [— (s, x)], where -d (s, x) is given by Eq. (174]) . 



APPENDIX F: EIGENFUNCTION WITH MINIMUM CHARACTERISTIC NUMBER AT s = 
At s = 0, Eq. (HIl) takes the form: 

To (x) = e-- 1 1 + ^/i (2^) e->i [To {y)\ | . (Fl) 

Let us differentiate both the parts of this relation with respect to x. It is easy to check the identity: 



d_ 

dx 



-^[e— ^/o(2V5^)]. 



After differentiating and using this identity, we integrate by parts on the right-hand side. The integrated term e ^ 
on the lower limit of integration, y — 0, will be cancelled by the result of differentiating the e^^. Finally, we have 

T^ (x) = e-^ r dylo {2^) e-y^\ [To (y)] T^ {y) . (F2) 
Jo 

Comparing this with Eq. at s = 0, we see that — Tq (x) is the (unnormalized) eigenfunction of this equation 

corresponding to the eigenvalue ^0 = 1- It is known from the theory of linear integral equations [25!| that the 
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eigenfunction corresponding to the maximum characteristic number can be chosen to be real and positive within the 
interval of integration. One can see that — Tq (x) < at any x > 0, and the corresponding eigenvalue /ig = 1 is indeed 
a maximal one. 
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